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Final  Report  Discussion:  Symmetry  Brealdng  Bifurcations  and 
the  Growth  of  Chaos  in  a  Rotating  Fluid 

ONR  Contract  N00014-86-K-(^85 
86  June  01  —  88  December  31 

Harry  L.  Swinney,  Principal  Investigator,  The  University  of  Texas  at  Austin 

Laboratory  experiments  and  numerical  simulations  have  been  conducted  on 
instabilities  in  flow  between  concentric  cylinders  (the  Couette-Taylor  system)  and 
Rossby  waves  and  vortices  in  eastward  and  westward  jets.  In  addition, 
Rayleigh'Benard  convection  near  the  onset  of  instability  has  been  investigated  in 
numerical  simulations  and  a  bifurcation  theory  analysis.  The  principal  results  are 
summarized  briefly  as  follows: 

Couette-Taylor  Flow  Instabilities  and  Tm*bulence 

1.  "Primary  instabilities  and  bicriticality  in  flow  between  counter-rotating 
cylinders,"  W.F.  Langford,  M.  Golubitsky,  R,  Tagg,  E.  Kostelich,  and  H.L.  Swinney, 
Phys.  Fluids  31, 776  -785  (1988). 

The  stability  of  the  basic  flow,  Couette  flow,  was  examined  in  a  linear  stability 
analysis.  The  analysis  method  was  not  new,  but  this  was  the  first  extensive  study  of 
the  primary  instability  for  a  wide  range  of  control  parameters  (cylinder  rotation 
rates  and  ratio  of  radii).  Particular  attention  was  focused  on  the  bicritical  curve., 
that  separate  the  transitions  from  Couette  flow  to  flows  with  different  azimuthal 
wavenumbers  m. 

2.  "Nonlinear  standing  waves  in  Couette-Taylor  flow,"  R,  Tagg,  S.  Edwards,  H.  L. 
Swiiuiey,  and  P.  S.  Marcus,  Phys.  Rev.  A39, 3734-3737  (1989). 

Recently  Demay  and  looss  conducted  a  nonlinear  stability  analysis  of  Couette 
flow  and  discovered  a  new  state  —  ribbons,  which  are  traveling  waves  in  the 
azimuthal  direction  but  standing  waves  in  the  axial  direction.  This  was  the  first 
new  primary  bifurcation  predicted  for  this  system  since  G.I.  Taylor’s  discovery  of 
axisymmetric  (Taylor)  vortices  and  spirals  in  1923.  We  have  observed  the  ribbons  in 
laboratory  experiments  and  corroborated  our  observations  with  direct  numerical 
simulations  at  the  parameter  values  where  the  ribbons  were  observed. 

3.  "Strange  attractors  in  weakly  turbulent  Couette-Taylor  flow,"  A.  Brandstater  and 
H.I,.  Swinney,  Phys.  Rev.  A  35, 2207-2220  (1987). 

The  transition  from  quasi-periodic  to  chaotic  Couette-Taylor  flow  was  studied  by 
laser  doppler  velocimetry  and  flow  visualization.  Velocity  data  were  used  to 
construct  attractors  and  Poincar^  sections.  The  attractor  dimension  was  found  to 
increase  from  2  at  the  onset  of  chaos  to  about  4  at  a  Reynolds  number  50%  above  the 
onset  of  chaos. 

4.  "Mass  transport  in  turbulent  Couette-Taylor  flow,"  W.Y.  Tam  and  H.L.  Swinney, 
Phys.  Rev.  A.  36, 1374-1381  (1987). 

Mass  transport  was  studied  at  very  high  Reynolds  numbers  R  (up  to  10®)  by 
injecting  a  pulse  of  dye  into  the  fluid  and  and  then  measuring  concentration  as  a 
function  of  time  and  axial  position.  Transport  in  the  axial  direction  was  found  to  be 
modeled  very  well  by  a  one-dimensional  diffusion  process,  with  the  diffusion 
coefficient  D  given  by  a  power  law,  D  R^,  where  (I  is  in  the  range  0.69  <  (i  <  0.86, 
depending  on  the  radius  ratio. 

5.  "Instabilities  and  chaos  in  rotating  fluids,"  II.  L.  Swinney,  in  Nonlinear 
Evolution  and  Chaotic  Phenomena,  ed.  by  G.  Gallavotti  and  P.  W.  Zweifel  (Plenum 
Publisliing  Co.,  1988),  pp.  319-326. 

This  brief  report,  contained  in  a  NATO  Advanced  Study  Institute  volume, 
summarizes  our  laboratory’s  recent  work  on  Couette-Taylor  and  geostrophic  flows. 


6 


Planetary-type  Flows 

6.  "A  laboratory  simulation  of  the  Great  Red  Spot  of  Jupiter,"  J.  Sommeria,  S.D. 
Meyers,  and  H.L.  Swinney,  Nature  331, 689-693  (1988). 

This  experiment  demonstrated  that  a  robust  large  coherent  vortex  can  form 
spontaneously  in  a  rapidly  rotating  turbulent  shear  flow.  The  fluid  was  contained 
in  a  rigid  rotating  annulus;  a  strong  westward  jet  was  generated  by  the  action  of  the 
Conolis  foi'ce  on  fluid  pumped  radially  outward.  The  vortices  that  form  in  this  flow 
exhibit  many  of  the  properties  of  jovian  spots:  they  have  the  same  sign  of  vortiHty  as 
the  background  flow,  they  move  along  with  the  surrounding  fluid,  the  width 
(longitude)  to  height  (latitude)  ratio  is  about  equal  to  two  for  the  large  spots,  and, 
most  importantly,  nearby  spots  are  attracted  to  one  another  and  merge.  Thus  this 
experiment  sheds  light  on  a  mystery  that  has  endured  for  more  than  three 
centuries,  the  origin  of  the  Great  Red  Spot. 

7.  'Numerical  simulation  of  Jupiter's  Great  Red  Spot,"  P.  S.  Marcus,  Nature  331, 
693-696  (1988). 

The  equations  governing  the  motion  of  an  inviscid  fluid  in  a  model  of  the  jovian 
atmosphere  are  solved  for  a  variety  of  initial  conditions.  Large  vortices  are  found  to 
form  spontaneously  in  chaotic  azimuthal  flows.  The  spots  are  stable  if  the  vorticity 
of  the  spots  has  the  same  sign  as  the  shear  of  the  surrounding  azimuthal  flow. 
These  simulated  spots  exhibit  many  of  the  properties  of  the  jovian  spots  and  of  the 
spots  observed  in  a  laboratory  experiment  (paper  number  6),  and  the  numerical 
simulation  yields  a  new  prediction  for  the  vertical  structure  of  the  Great  Red  Spot. 

8.  "Laboratory  model  of  a  planetary  eastward  jet,"  J.  Sommeria,  S.D.  Meyers,  and 
ILL.  Swimiey,  Nature  337, 58-61  (1989) 

A  strong  eastward  jet  is  generated  in  the  rotating  annulus  (used  in  paper 
number  6)  by  reversing  the  direction  of  the  radial  flow.  For  a  wide  range  of 
parameters  Rossby  waves  form  on  a  narrow,  essentially  laminar,  eastward  jet.  The 
wave  velocity  and  wavelength  are  accounted  for  by  simple  scaling  arguments. 
Experiments  with  dye  injected  into  the  annulus  far  from  the  jet  show  that  the  dye  is 
rapidly  mixed  away  from  the  core  of  the  jet,  but  the  jet  acts  as  a  strong  barrier  to 
tracer  transport;  photographs  obtained  hundreds  of  rotations  after  the  injection 
dramatically  illustrate  the  sharp  barrier  formed  by  the  jet  core.  This  could  have 
important  implications  for  the  understanding  of  transport  of  pollutants  in  the 
atmosphere  and  oceans. 


Convection 

9.  "Global  bifurcation  to  traveling  waves  in  axisymmetric  convection,"  L. 
T^ckerman  and  D.  Barkley,  Phys.  Rev.  Lett.  61, 408-411  (1988). 

We  have  discovered  a  traveling  wave  state  in  axisymmetric  convection  in  which 
rolls  are  continually  created  at  the  sidewalls  and  annihilated  at  the  center. 
Heuristically,  for  an  aspect  ratio  of  five,  the  system  alternates  between  four  and  five 
rolls  in  an  irreconcilable  wavelength  selection  conflict.  We  have  demonstrated  that 
for  conducting  side  walls  this  oscillatory  state  is  a  heteroclinic  orbit:  a  saddle-node 
bifurcation  that  leaves  in  its  wake  a  limit  cycle  whose  period  is  infinite  at  onset.  In 
contrast,  for  insulating  sidewalls,  when  the  central  roll  disappears,  again  via 
saddle-node  bifurcation,  transition  occurs  instead  to  a  steady  four-roll  state.  We 
have  interpreted  the  Prandtl  number  dependence  of  the  transition  threshold  in  light 
of  a  proposed  wavelength  selection  criterion  [Y.  Pomeau  and  P.  Manneville,  J. 
Phys,  (Paris)  42, 1067  (1980)]  prohibiting  net  radial  flow  in  axisymmetric  convection. 
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1 0.  'Traveling  waves  in  axisymmetric  convection:  the  role  of  sidewall  conductivity," 
D.  Barkley  and  L.  Tuckerman,  Physica  D,  in  press  (1989);  see  also  the  paper  in  New 
Trends  in  Nonlinear  Dynamics  and  Pattern  Formation:  the  Geometry  of 
Nonequilibrium,  NATO  Advanced  Study  Institute  Series,  ed.  by  P.  Cuullet  and  P. 
Huerre  (Plenum,  New  York,  1989). 

Traveling  waves  in  axisymmetric  convection  such  as  we  have  computed 
numerically  have  been  predicted  by  theory  to  exist,  but  have  never  been  observed 
experimentally.  We  have  explained  this  by  showing  that  the  four-roll  branch, 
which  intercepts  the  traveling  wave  solution,  is  an  isola  whose  range  of  existence 
increases  as  the  ratio  of  sidewall  to  fluid  conductivity  is  decreased.  Traveling  waves 
will  be  triggered  to  occur  only  if  the  conductivity  of  the  sidewalls  is  over  60  times  that 
of  the  fluid. 
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and  U.  Schumann,  in  Proceedings  of  the  Third  GAMM  Conference  on  Numerical 
Methods  in  Fluid  Mechanics,  ed.  by  E.  H.  Hirschel  (Vieweg  and  Sohn, 
Braunschweig,  1980)]  as  algebraic  equations  relating  coefficients  in  one  nonperiodic 
direction.  We  have  developed  a  matrix  formalism  for  a  general  geometry,  in  which 
the  influence  matrix  method  can  be  recast  as  an  application  of  the  classic 
Sherman-Morrison-Woodbury  formula  of  numerical  linear  algebra.  This 
formalism  has  permitted  us  to  apply  the  method  to  cylindrical  coordinates, 
complicated  by  boundaries  in  two  directions  and  by  the  coordinate  singularity  on  the 
axis. 
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Implicit  time-stepping  in  spectral  methods  required  the  inversion  of  elliptic 
operators  acting  on  orthogonal  function  expansions.  This  operation  is  generally 
quadratic  in  the  number  of  basis  functions  in  a  nonperiodic  (bounded)  direction,  but 
can  be  made  linear  in  at  most  one  Cartesian  direction  by  use  of  a  well-known 
recursion  relation  [D.  Gottlieb  and  S.  A.  Orszag,  Numerical  Analysis  of  Spectral 
Methods:  Theory  and  Applications  (SIAM  Press,  Philadelphia,!  977)].  We  have 
isolated  the  relevant  property  responsible  for  the  existence  of  such  recursion 
relations,  and  applied  it  to  spherical  and  cylindrical  geometries. 
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Steady-state  continuation  via  Newton's  method  is  an  attractive  alternative  to 
time  stepping,  both  because  dynamical  timescales  need  not  be  resolved,  and  because 
unstable  steady  states  are  accessible.  However,  Newton's  method  requires  the 
solution  of  large  linear  systems  involving  Jacobian  matrices  which  are  full, 
nonsymmetric,  and  poorly  conditioned.  We  describe  a  way  of  modifying  a 
time-dependent  code  to  perform  Newton  iteration,  in  which  implicit  time  stepping 
automatically  acts  as  an  effective  preconditioner  for  the  Jacobian  matrices.  We  have 
been  able  to  compute  steady  states  in  l/30th  the  time  required  by  a 
Courant-condition-limited  time-stepping  code  by  using  the  NSPCG  biconjugate 
gradient  squared  algorithm  [T.  C.  Oppe,  W.  D.  Joubert,  and  D.  R.  Kincaid  (1988), 
Publ.  3216  of  the  Center  for  Numerical  Analysis,  University  of  Texas  at  Austin]. 


